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Abstract 

Phase-locked solutions of coupled oscillators are studied with asym- 
metric coupling strengths or inhomogeneous natural frequencies. The so- 
lutions show remarkable profiles of phase lags from the pacemaker corre- 
sponding to the ratio of upward and downward coupling strengths. By 
means of the existence condition of phase-locked solutions, the transition 
points from phase-locked to phase drift states are estimated. The applica- 
tion of the existence condition to the case of the linear gradient of natural 
frequency illustrates some scaling properties in the frequency diagrams. 

1 Introduction 

Coupled oscillators have been dealt with in wide and diverse areas associated 
with oscillation and synchronization phenomena. The phenomenon of collective 
synchronization of phases or frequencies attracts attentions in not only biology 
but also in physics and engineering fj],0]. Theoretical studies have been deeply 
carried out for the case of global coupling, and clarified characters of phase 
transition of entrained oscillators [prQ- The case of local coupling has also 
been studied in various contexts motivated by rhythmic phenomena in biological 
systems. They include animal gaits B, fish swimming [P,|lC|], and the peristaltic 
movement of gastrointestinal tracts |ffl, [ll], [l^] . In recent years, researches of 
coupled oscillators are extended to the case of non-local (but not global) coupling 
]l3| and more complicated network systems Q. 

Since such complicated networks have, in general, intricate boundaries and 
inhomogeneous couplings, we have been found little clues to analyses of dynam- 
ical behaviour of them. We can, however, analyse some phase-locked solutions 
for more simple and symmetric trees of couplings. In the case of m-ary trees 
with the height n and to™ leaves, for instance, the system of coupled oscillators 
has some particular solutions with identical phases in each level, and these so- 
lutions are reduced to the phase-locked solutions for chains of oscillators with 
asymmetric coupling strengths Eq] . As well as complicated network structures, 
intricate dynamical processes appear in biological systems, of which it is impos- 
sible to drag the inherent or isolated parts out. It is one of the ways to study 
such systems that we construct models from the phenomenological evidence and 
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draw qualitative features of them. When we compare the models with data ob- 
served in real systems such as efficiency of movement and the transport volume, 
we will need to know the dynamical behaviour and solutions of the models in 
detail. 

In this paper, we study chains of coupled phase oscillators. It is the remark- 
able point that we can derive phase-locked solutions of the systems by algebraic 
calculations, even if the coupling strengths are asymmetric or the natural fre- 
quencies are inhomogeneous. Moreover, the analytical form of the solutions is 
not much cumbersome, and thus the profile of solutions are easy to see. From 
the existence conditions of the phase-locked solutions, we can estimate transition 
points of phase-locked to phase drift states. 

We present coupled phase oscillators with asymmetric coupling strengths in 
the following section, and derive the phase-locked solutions in section 0. The 
profile of the phase-locked solutions and bifurcation phenomena are studied 
for the cases of asymmetric coupling strengths (section ^) and of a linear fre- 
quency gradient (section ||). In section |6[ we discuss the difference analogue of 
the continuous system corresponding to the oscillator systems with asymmetric 
coupling strengths. 

2 Coupled phase oscillators 

A linear chain of phase oscillators with asymmetric couplings is governed by the 
following equations of phases 9j (J = 0, . . . , n): 

eo = wo + a*h(6i-B a ), 

9j = ujj + a d h(0j-i - 9j) + a u h(9 J+1 - 9j), (j = 1, . . . , n - 1), (1) 
8 n = w n + a* d h(9 n _i - n ), 

where u>j is the natural frequency, a u and ad are coupling coefficients, and h is a 
periodic function. For the Oth and nth oscillators, coupling coefficients, a* and 
a* d are set according to the boundary conditions. 

Introducing new variables defined by neighbouring sites, 

ipj := 9j - 9j-i, dj := lu 3 - uij-i, 

we obtain n equations for the phase differences: 

ip = d + A u h+-A d h-, (2) 
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In the analysis for phase locked/drift solutions shown below, we assume h is 
an odd function, h(—ip) = —h(ip), and thereby h = h + . 



3 Phase lags 

Phase-locked solutions of eq. |(2)| mean all the oscillators have the same period 
and the phase difference is constant, ipj =0. Thus they are solutions of 



Ah + d = o, 



(3) 



where A := A u + Ad and h := h + . We solve formal solution of the above 
equations under the free boundary conditions (mutual cntrainment), and one- 
sided coupling conditions (forced entrainmcnt). 

Fast, we study the case of free boundary conditions, a* = a u , a* d = ad, where 
the boundary oscillators and their neighbours are coupled in the bi-directional 
manner. Under these conditions, all of the oscillators have the same period 



9j — Sl n / A r 



(4) 



where fl n and A„ are defined below. From eq. (3), we get the solution hj 
(j = 1, . . . , n), the jth element of h, in the form of 



a d hj 



(A„/A n )A„_j 



or 



where 



A 
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(A n /A n )A 



(5) 
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Aj := y^(u n -j +k - uq) ( — ) , Aj := ^(w„ - Wj.fe)^ — 
^ ad ' . \ a*. 



k=0 



k=0 



The phase difference, ipj, is obtained by inverting eq. |(5)| . When the coupling 
function h is continuous and bounded, solutions exist only if min h < hj < max h 
for all j. The linear-stability conditions of solutions are obtained from the 
eigenvalue problem of the linearized matrix, Ah' , where the jth elements of h' 
is h'(?pj)- 

We comment the cases of forced oscillations, that is, one of the terminals is 
a forced oscillator. Suppose that the 0th oscillator is one-sided coupling to its 
neighbour, a* = 0, but the opposite end satisfies the free boundary condition, 
a* d = a d . Under these conditions, a solution of eq. (3) is hj — A n -j/ad- If 



the boundary conditions are exchanged, that is, the nth oscillator is a forced 
oscillator, a* d = 0, and the 0th oscillator satisfies the free condition, a* = a u , a 
solution of eq. (3) is obtained in the form of hj — Aj_i/a u . 
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4 Asymmetry couplings 



We study the mutual entrainment to a pacemaker with asymmetry couplings in 
the present section. We assume that all the oscillators have the same natural 
frequency but only the Oth oscillator, the pacemaker, has the high frequency If 
we set uj Q > and u>x = ■ ■ ■ = u n = 0, eq. 



5j becomes 



The entrained frequency is obtained from eq. (4) 



f w /(l + --- + A"), (a u >a d ), 
w n = < (8) 
\ w A n /(l H h A"), (a u <a d ), 

where A := m.ax(a u ,ad)/ min(a u ,ad). We set a u , ad, and h'(0) are positive be- 
low. For numerical calculations shown in figures, we chose the coupling function 
as h(tp) = s'mip + 0.1 sin 2^ — 0.03 sin 3ip to obtain generic results. 

When the coupling strengths are asymmetric, a u ^ ad, the profile of hj 
given in eq. |7J has the exponential dependence on j. In fig. |l|, we illustrate 
some features of phase- locked solutions: the profile of coupling interaction, phase 
difference and phase lag from the pacemaker. Calculations of these quantities 
are carried out for the system of 20 oscillators (n = 19) and for three ratios of 
couplings, a u /ad = 2, 1, and 1/2. If the coupling function h(ip) is approximated 
to the linear function for small ip, then the dependence of the phase lag, 9j — 6q, 
on j is almost (i) exponential for a u > ad, (ii) parabolic for a u — ad, and (hi) 
linear for a u < ad- Such dependence is inherited for more complicated network 
systems on trees fl5|| . 

For small coupling strengths, the phase- locked state is broken and the phase 
of each oscillator drifts from the phase of the pacemaker. Such a phase drift 
solution has two regions of the frequency-entrainment since the pacemaker and 
other oscillators have different natural frequencies. The phase drift occurs be- 
tween these two regions. Fig. ^| shows the transition between phase-locked and 
phase drift solutions for the three ratios of a u and ad- Each system has 10 os- 
cillators (n = 9). The averaged frequency of the jth oscillator and the coupling- 
strength are defined by 

(uij) := lim (0j(t) — 0j(O))/t, e := min(a u , ad) max h. (9) 

t— >oo 



In fig. ||, (ujj) is obtained from numerical integrations of eq. [1] with the step 



width in e, Ae = 10~ 3 x max/i. We can estimate the transition point of phase- 
locked and phase drift solutions from the necessary condition for the existence 
of the phase-locked solutions, 

1 + • • • + A"- 1 
e > e n := u - 



! + •■• + A' r 



The critical point (e„,d>„) is plotted by the open circle in fig. 0. 

Numerical calculations show that phase drift states always consist of two re- 
gions: the pacemaker and the other oscillators. It implies that if the pacemaker 
can entrain its neighbouring oscillator, all oscillators will be entrained to the 
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pacemaker, but if not so, no oscillators entrained to the pacemaker. It is also 
possible to estimate this entrainment feature from the existence condition of the 
phase-locked solutions. We assume that the pacemaker is entraining j oscilla- 
tors. Then we divide the system into two sets of coupling oscillators: (i) upper 
(j + 1) oscillators including the pacemaker, and (ii) lower (n — j) oscillators. For 
the upper set (i), the phase-locked solution exist only if 

For the lower set (ii), we add one oscillator with the entrainment frequency of 
the set (i) to the upper side as the pacemaker. If the coupling strength e satisfies 

_ l + '-' + A™-^ 1 

e<uJ-j— ■ — , (11) 

3 H h A"-J 

then the phase drift solution exists for the set (ii) with the pacemaker. Since 
no value of the coupling strength e satisfies both of conditions |(10)| and (11) 
oscillators entraining to the pacemaker are all or none. 



5 Linear frequency gradient 



Daido 1 12 1 has studied chains of coupled phase oscillators with a linear gradient 
of natural frequencies by numerical calculations. He has shown some scaling 
properties in the frequency diagram, a plot of averaged frequencies against the 
coupling strength. In the present section, we attempt to derive some of them 



from eq. |[5) 



We assume that the natural frequency of the jth oscillator is u)j := j/n and 
coupling strengths are symmetric a := a u = ad. Then a phase-locked solution 



is obtained from eq. [5] in the form of 

hj = 6 -^±j^l, (12) 

where S :— l/n is the difference of natural frequency between neibouring oscil- 
lators. The entrained frequency of this solution is 

w, =5—1 = 



We set a and h'(0) are positive, and chose the coupling function as the previous 
section for numerical calculations shown in figures below. 

We note that hj takes the maximal value at the centre of the chain of en- 



trained oscillators from eq. (12). If the phase difference ipj is small enough to 
approximate that the coupling function h(ip) has the linear dependence on ip, 
ijjj has the convex dependence on j as well as hj. Then the phase lag (6j — 6q) 
has an inflection point at the centre of the chain. Figure [| shows the profile of 
phase lags and other variables. 

When the coupling strength is diminished, the entrained region is divided 
into some clusters of frequency-entrained oscillators, called frequency plateaus 
[Q,|l^,|l2|. In the limiting case that the coupling strength is vanishing, entrained 
domains of oscillators will split into each oscillator. 
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Although the frequency diagram shows fine and complicated bifurcation 
structure |l2] ] , we can estimate the approximate arrangement of transition points 
by means of the necessary conditions for existence of phase- locked solutions. To 
illustrate the frequency diagram, we define the averaged frequency of the jth os- 
cillator and the coupling strength as eq. |(9)| . We pick out the frequency plateau 
with m oscillators from the (n + l)-oscillator system, and consider these oscil- 
lators as the isolated system with the free boundary conditions in both sides. 
Phase-locked solutions of the m-oscillator system exist only if 

{(S/8)m 2 , (to : even), 

(14) 
(S/8)(m 2 - 1), (m:odd). 

The entrained frequency is the average of natural frequencies, 

&k,k+m-l = (Wfe + • • • + Ldk+m-l)/™ 

= 5(m + 2k- l)/2, (fe = 0,...,n-m). (15) 

In fig. ^, we put the approximate critical point (ujk,k+m-i, £ m ) on the fre- 
quency diagrams for the system of 10 oscillators (n = 9). Some of approximate 
points are located near the transition points, but others have no transition points 
around them. The latter points indicate that it is impossible for the system to 
have some states such as small plateaus located near the boundaries or ex- 
tremely asymmetric arrangements of plateaus. We obtain the scaling behaviour 
of 

^0,771—1 

to e m /e n+ i from eqs. pi] and |(I5)1 as 



1/2 

Wo.m-l ~ ^ ( — — ) (16) 

1 \ e n+i. 



in the limit of n — > oo. 



6 Discussion 



We can consider the phase-difference system (2) is derived from some contin- 
uous partial differential equation by means of the difference analogue. If the 
coupling function h is assumed to be the odd function, the continuous system 
corresponding to eq. (2) becomes 

d t ijj + Td x h{^) = d{x)+Ddlh{iP), (17) 

where x is the spatial coordinate. Coefficients T and D are defined by 



T := (a<j - a u )Ax, D 



— (Ax)" 2 



for the difference interval Ax. Here we took the central difference for the dif- 
ferential of first order. In the continuous system (17), entrained solutions are 
obtained from the ordinary differential equation in the form of 



d 2 h dh „ . 
^- r d^=°- 



(18) 
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When all of natural frequencies are identical, d(x) = 0, solutions of eq. |(18)| 
J c+bexp(Tx/D), (T ^ 0), 

[ c + bx, (r = o), 



are 



where b and c are co nstan ts defined by boundary conditions. These solutions 
are consistent with eq. (7] and fig. [l]. In another case co rresp onding to the linear 
frequency gradient, d(x) = S (const.), a solution of eq. (18) are obtained in the 
form of 

h(x) =c + bx- {6/2D)x 2 . 



This is also consistent with eq. (12) and fig. |3|, if we set arbitrary constants b 
and c in proper values. 

Here we give more examples. When the both cases mentioned above are 
combined, that is, the natural frequency has the linear gradient in space, d{x) 
S, and coupling strengths are asymmetric, f^O, a solution of eq. |(1£ 



is 



h(x) =c+ (TS/D)x + bexp{Tx/D). 

Even if r = 0, we can obtain the similar dependence of h on x from the system 
in which the natural frequency has the exponential gradient in space as d(x) = 
7exp(7(x — L)) (L is the system size). A solution of this system is 

h(x) = c + bx — (I/D7) exp(7(x — L)) . 

In the last, we comment on the case that the coupling functions is divided 
into the odd and even parts. In general, the division of a function h{x) is carried 
out in the form of 

h(x) = h odd (x) + h eve n(x), 

_ h(x) - h(-x) _ h(x)+h(-x) 

^odd^J • 1 ^cvcn\X ) . 



The system of coupled phase oscillators |(2)| is derived from the following con- 
tinuous system by the discretization of space 

d t ip + r odd a K /i dd(V') - 1 even d x 1 ''even 

= d(x) + D odd dlh od d(ip) - D cvcn dlh CVC n(ip), (19) 

where coefficients are defined as 

Todd := {ad - a u )Ax, T cvcn := (ad + a u )Ax 

a d + a u , , 2 a d -a u 2 
-Dodd := 2 ' Dcvcn := 2 ^ ' ' 

We note that the coefficients of advection and diffusion terms are almost ex- 
changed between the odd and even parts, 



Axid = (Aa;/2)r even , = (Ax/2)r o 



dd- 



It is also possible to obtain the phase-locked solutions from eq. |(19)| as above, 
but we do not mention them here. 
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7 Conclusions 



We obtained the phase-locked solutions of coupled phase oscillators with asym- 
metric coupling strengths or inhomogeneous natural frequencies. These solu- 
tions show the different profiles of phase lags between the three types of ratios 
of the upward to downward coupling strength. The phase differences obtained 
in analytical form, must be applicable to consider the efficiencies of biological 
motions and transportations based on the collective oscillation. Moreover, in- 
herent structures of the real systems, the asymmetry of couplings, for example, 
are suggested by the spatio-temporal patterns observed in the systems. 

From the existence conditions of phase-locked solutions, we derived transi- 
tion points from phase- locked to phase drift states. We also obtained that the 
length of entrained region to the pacemaker is approximately zero or the sys- 
tem size. However, for the discrete time-dependent Ginzburg-Landau equations, 
the number of entraining oscillators to the pacemaker depends on the coupling 
strength, and the sequence of bifurcations is observed numerically Jl5]] . 

Next, we showed the profiles of phase-locked solutions for the linear fre- 
quency gradient. The phase differences take the maximal absolute value around 
the centre of the chain system. This means that the phase drift tends to oc- 
cur around the centre when the coupling strengths are diminished. As well 
as the case of the asymmetric coupling strengths, we estimate the approxi- 
mate transition points in the frequency diagram from the existence conditions 
of phase-locked solutions, and derive some scaling properties. 
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Figure 1: The profiles of phase-locked solutions for phase oscillators with 
asymmetric couplings: the coupling interaction, hj, phase difference, ipj, and 
phase lag from the pacemaker, (8j — do), where ipj and (8j — 8q) are plotted 
with 7T units. The left, centre, and right columns are the cases of (i) a u /a d = 2 
(a u = 1.0, a d = 0.5), (ii) a u /a d = 1 (a„ = 1.5, a d = 1.5), and (iii) a u /a d = 1/2 
(fflu = 1.5, a d = 3.0), respectively. Solid lines are connecting the points derived 
eq. |(7)| , while open circles are obtained from numerical calculations of the phase 
equation |(1)| . 
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Figure 2: The frequency diagrams for phase oscillators with the asymmetric 
couplings. The left, middle, and right diagrams are the cases of (i) a u /a,d = 2, 
(ii) a u /a e i = 1, and (hi) a u /a e i = 1/2, respectively. The averaged frequencies 
(ujj) are plotted against the coupling strength e. The dots are obtained from 



numerically calculations of eq. (1). The open circles denote critical points of 



the existence condition for phase-locked solutions. 



l.Ofl 
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Figure 3: The profiles of phase- locked solutions for phase oscillators with linear 
frequency gradient: the coupling interaction, hj, phase difference, -0^, and phase 
lag from the Oth oscillator, (6j —Oq). The upper left shows the natural frequency 
of each oscillator. 
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Figure 4: The frequency diagram for phase oscillators with the linear gradient 
of natural frequencies. The diagram is plotted the averaged frequencies (uj) 
against the coupling strength e. Averaged frequencies plotted by dots are nu- 
merically calculated. The open circles denote critical points for the existence 
condition of phase- locked solutions with m oscillators (m = 1, . . . , n + 1). 
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